Two-color discrete localized modes and resonant scattering 
in arrays of nonlinear quadratic optical waveguides 

Mario I. Molina 1 , Rodrigo A. Vicencio 1,2 , and Yuri S. Kivshar 3 
1 Departamento de Fisica, Facultad de Ciencias, 
Universidad de Chile, Casilla 653, Santiago, Chile 
2 Max- Planck- Insitut fur Physik komplexer Systeme, 
Nothnitzerstr. 38, Dresden 01187, Germany 
3 Nonlinear Physics Center, Research School of Physical Sciences and Engineering, 
Australian National University, Canberra ACT 0200, Australia 

Abstract 

We analyze the properties and stability of two-color discrete localized modes in arrays of channel 
waveguides where tunable quadratic nonlinearity is introduced as a nonlinear defect by periodic 
poling of a single waveguide in the array. We show that, depending on the value of the phase 
mismatch and the input power, such two-color defect modes can be realized in three different 
localized states. We also study resonant light scattering in the arrays with the defect waveguide. 
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I. INTRODUCTION 



The study of nonlinear propagation of light in periodic photonic structures recently at- 
tracted strong interest due to the unique possibility of observing experimentally an interplay 
between the effects of nonlinearity and periodicity, including the generation of spatial opti- 
cal solitons Q|. Recently fabricated nonlinear periodic structures such as arrays of weakly 
coupled nonlinear optical waveguides can simultaneously support distinctive types of self- 
trapped optical beams because linear effects such as diffraction may differ dramatically 
compared to those in the corresponding continuous systems [2-5]. 

In addition to the nonlinear periodic systems, during last years a growing interest is ob- 
served in the study of nonlinear optics associated with the so-called quadratic nonlinearities 
which may produce the effects resembling those known to occur in cubic nonlinear mate- 
rials. Typical examples are all-optical switching phenomena in interferometric or coupler 
configurations as well as the formation of spatial and temporal solitons in planar waveg- 
uides 0. One of the recent highlights in this field is the first experimental demonstration of 
discrete solitons with two frequency components mutually locked by a quadratic nonlinear- 
ity 0. These optical experiments have been performed in arrays of weakly coupled channel 
waveguides with tunable cascaded quadratic nonlinearity, and they demonstrated a good 
agreement with the theoretical analysis 3- As a matter of fact, it' was demonstrated that 
arrays of coupled channel waveguides fabricated in a periodically poled Lithium Niobate slab 
represent a convenient system to verify experimentally many theoretical predictions. These 
experimental observations open many novel perspectives for employing larger nonlinearities 
in the lattice systems with quadratic materials. 

In this paper, we analyze the properties and stability of two-color discrete localized modes 
in arrays of channel waveguides. In particular, we show that when periodic poling is applied 
to just a single waveguide in the array, it creates a nonlinear defect 0, Q] that may support 
strongly localized discrete modes similar to discrete solitons and also display specific resonant 
scattering properties. 

The paper is organized as follows. In Sec. II we introduce our model and find the profiles 
of two-color discrete localized modes. We show that such modes can exist in three different 
states, and they can be observed in the same array depending on the array parameters 
and the value of the input power. Section III is devoted to the study of the resonant light 
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transmission through the defect waveguide that can be associated with the Fano resonance. 
Finally, Sec. IV concludes the paper. 



II. TWO-COLOR LOCALIZED MODES 

Being motivated by the design of the periodic photonic structures recently employed 
for experiments [jj, we consider an array of weakly coupled linear waveguides where one 
waveguide has periodic poling, and therefore it possesses a quadratic nonlinear response. 
When the matching conditions are satisfied, the fundamental-frequency (FF) mode with 
the frequency u generates the second-harmonic (SH) wave at the frequency 2ou, so that 
such a structure with the poled waveguides may behave as a nonlinear defect with localized 
quadratic nonlinearity. 

In the tight-binding approximation the effective equations for the complex envelopes 
of the FF wave (u n ) and its SH component (v n ) coupled at the defect waveguide (at n — 0) 
can be written in the form 

" + c u (u n+1 + u n -i) + 2ulv 5 n0 = 0, 



dv 

i—j^- + c v (v n+ i + u n _i) - Av n + u 2 5 n Q = 0, 

where c u and c v are the coupling coefficients, A is the phase mismatch parameter, and S n Q 
is a delta-function. 

The discrete model ((TJ) has two conserved quantities, the total power 

P= (l^| 2 + 2|t; n | 2 ) (2) 

n=— oo 

and the Hamiltonian. We look for stationary solutions of the coupled equations in the form 

u n {z) = u n e iXz , v n (z) = v n e 2iXz , 

where we assume, without loss of generality, that the amplitudes u n and v n are real. As a 
result, Eqs. (JT} become 

Xu n = c u (u n+1 + u n _i) + 2u v 5 n0 , 

(3) 

2Xv n = c v (v n+1 + w n _i) - A v n + u 2 Q S n0 . 
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We look for spatially localized discrete states in the form u n = Uq 7 ? an d v n = t>o£' n ', 
where \r]\, |£| < 1, and obtain from Eqs. © at n ^ two relations, A = c u (r] + r/ _1 ) and 
2A = c„(£ + — A, which can be combined to yield 

c v U + 7 J - 2c„ h + - J = A- (4) 



On the other hand, for the defect waveguide at n = we obtain A = 2c u r] + 2u and 
2At>o = 2c„t>o£ _ At; + it§. After some simple algebra, we solve for the parameters uo and 
t>o in terms of rj and £ and obtain 

= y Q - rjj , = c^ Q - £J ■ (5) 
In addition, Eq. (J2J) can be expressed as 

and the next step is to replace the amplitudes uq and vq in Eq. © and obtain the second 
equation, 

Equations (J3| and (0) constitute a system of two coupled nonlinear equations for the func- 
tions r\ and £. Moreover, Eq. depends on the mismatch parameter A, whereas Eq. (J7|l 
depends on the total power P. 

First, we assume that c u , c v > 0. From Eq. (jSJ) we find that if v® > then < £, r/ < 1. 
On the other hand, if vq < 0, then — 1 < £,77 < 0. Also, without loss of generality, we 
can consider only Uq > 0. Thus, we should find real solutions for the values 77 and £ in the 
domain defined by the conditions < 77, £ < 1 and — 1 < rj, £ < on the plane (77, £). 

From the structure of Eqs. (jlj) and (Jlj) we notice that a change A — > —A is equivalent to 
the change (77, £) — > (—77, — £). Thus, we only need to consider positive values of A. In order 
to find spatially localized states, we should analyze the functions of 77 and £ defined by the 
left-hand sides of the expressions (j3J) and (|7|1. for varying values of the mismatch parameter 
A and the total power P. This analysis reveals that, in general, for a given value of A, a 
minimum value of the power is required to support a localized state at the defect. 

As an example, in Figs, l(a-d) we show several profiles of the localized modes for the 
case A = and two values of the power P, for the typical values c u = 1 and c v = 0.5. 
Below a certain threshold power, there exist no localized states, but above the threshold 
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FIG. 1: Examples of localized modes at A = and two values of the total power above the 
threshold: (a,b) P = 5 and (c,d) P = 10. Solid and dashed curves show the fundamental-frequency 
(u n ) and second-harmonic (v n ) fields , respectively. 





FIG. 2: Same as in Fig. 1 but for A = 5 and two values of the total power above the threshold: 
(a,b) P = 7 and (c-e) P = 10. In this case, two staggered localized modes appear about the 
threshold power, and then additional unstaggered mode appears for larger values of P. 

two localized modes appear, unstaggered (a) and staggered (b) modes, respectively [see 
Figs. l(a,b)]. As the total power P increases, these two modes become more localized, as 
shown in Figs. l(c,d). 

For a finite positive mismatch A and increasing power P, we find that there is a threshold 
power below which no localized states exist. When the power exceeds the first threshold, at 
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least one (staggered) localized modes becomes possible and then the other mode appears for 
larger powers. As the value of the mismatch A increases further, it can be proven that the 
contour curve associated with Eq. (J3J), which depends on A, experiences a curvature change 
at A = 2(2c u — c v ) = 3 and, as a result, a double root is possible at a certain power level. 
Further power increase gives rise to two staggered localized states. Even further increase in 
the power makes possible another (unstaggered) localized state to exist, as shown in Fig. 2 
at A = 5. 

The threshold value of the total power that corresponds to the appearance of the first 
localized mode can be found analytically in the form, 

P mm (A)=c u [(A + 4cJ 2 -4a 1 / 2 , (8) 

for the condition A > 2(c v — 2c u ), and 

^L(A) = c4(A-c u ) 2 -4a 1 / 2 , (9) 

for the condition A < 2(2c u — c v ). In the first case, the lowest localized mode is unstaggered, 
whereas in the second case, it is staggered. Applying these results to Fig. 1, we obtain 
^ mi „(0) = J P min (0) = v / 15 = 3.873. 

By solving Eq. ((TJ) numerically, we have checked the stability of all localized modes to 
propagation, and we came to the conclusion that such modes are stable, and they do not 
display any bistability. 



III. RESONANT TRANSMISSION 



The two-color localized modes described above can be excited and generated experimen- 
tally in arrays of weakly coupled quadratic waveguides being detected through the specific 
features of the transmission coefficient Q]. Therefore, here we analyze the scattering of a 
plane FF wave by a quadratic defect waveguide. When the phase-matching conditions are 
satisfied, after the interaction with the quadratic waveguide, the FF wave generates a SH 
wave which could either propagate or get trapped being guided by the defect waveguide in 
the form of a localized mode. 
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Propagating SH field. To calculate the transmission coefficient in the case of the propa- 
gating SH field, we present the fields as 

ikn _|_ e ~ikn < q 

u n (z) = e^ z { (10) 
cie ifcn , n > 0, 

Vn(*)=e^ (11) 
c 2 e iqn , n > 0, 

where a is the amplitude of the FF wave before the scattering, and c\ )2 are the amplitudes 
of the reflected and transmitted FF and SH waves, respectively, and the wave numbers k 
and q are defined in the domain < k, q < it. Using these expressions far from the defect 
site n = 0, from the phase-matching condition X v = 2X U we obtain 

A . 
c v cos q — — = 2c M cos k. (12) 

Equation (fT2"j) defines q = q(k) which is real and positive in some interval, k min < k < /c max , 
where 

-i f A±2c v 

Outside this interval, q is purely imaginary, and this case corresponds to the SH field localized 
at the waveguide. 

Evaluating the fields at the sites n = —1,0, we obtain 

A n ci = c u {cie ik + ae~ ik + he ik ) + 2c\c 2 , 

(13) 

Kc 2 = c v (c 2 e ig + b 2 e iq ) - c 2 A + c\, 

Ci = a + &!, b 2 = c 2 , (14) 

the latter condition implies that the SH fields are generated symmetrically. Using these 
results, we find the important relation c 2 = (i/2c v sin q)c\ that allows us to obtain the non- 
linear equation for the transmission coefficient of the FF wave, defined as t(k) = |ci| 2 /|a| 2 , 

m = [i + B^mr (15) 

where B(k) = |a| 2 A(/c) = \a\ 2 /{2c u c v sin k sinq(k)}. By rewriting Eq. (|15jl in the form 

A 2 (k)t 3 + 2B{k)t 2 + t = 1, (16) 



and using the fact that A(k) > 0, it is easy to see that this result does not predict bistability. 
As a matter of fact, the only real and positive solution of Eq. (fT3j) can be found in an 
analytical form from the equation 

Localized SH field. To calculate the transmission coefficient in the case of the localized 
SH field, we present the SH field as 

iXx (b 2 e 9n , n<0, 
v n = e tXvZ { (17) 
I c 2 e- qn , n>0, 

which can be obtained by replacing q — > iq in Eqs. (fTTj) . Real and positive values of q can 
be found from Eq. (|T2|) that now takes the form: 

c v cosh q — — = 2c u cos k. (18) 

Wavenumber k of the incident FF field should be taken outside of the domain (/c m i n , k max ), 
i.e. for < k < k m i n and /c max < k < it. The relation between the transmitted amplitudes 
now becomes c 2 = c\/(2c v sinh q), and the transmission coefficient t(k) is defined by the 
equation 

m = i + z)WW (19) 

where D{k) = |a| 2 /[2c„c u sin A;sinh(j»(fc)]. Equation (fTHjl possesses areal and positive solution 
for t(k) that can be found in an analytic form. 

Figures 3(a-d) show some examples of the dependence of the transmission coefficient t(k) 
on the wavenumber k (0 < k < tt) for different values of the coupling coefficients c„, c v 
and the mismatch parameter A. In a sharp contrast with the problem of the cubic defect 
analyzed earlier in Refs. where the transmission coefficient does not show any interesting 
structure and possesses a single maximum only, here the transmission coefficient shows two 
minima associated with the resonant suppression of the transmission. For small values of 
the ratio c v /c u , the middle interval is narrow and the transmission remain low. The effect of 
the nonzero mismatch parameter [see Fig. 3(d)] is to shift the position of the middle interval 
making the transmission curve asymmetric. 

The resonant suppression of transmission at some points, i.e. t(k m i n ) = t(k max ) = 0, 
corresponds to a novel type of the well-known Fano resonance Q| as recently discussed in 



Ref. 



12| for a simpler model. Indeed, destructive interference and resonant suppression of 



transmission is observed when there exists an extra localized state coupled to the propagating 
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FIG. 3: Transmission coefficient of a plane wave (shown for the fundamental frequency only) after 
scattering by a nonlinear quadratic defect, for A = and (a) c v = 1, c u = 0.1, (b) c v = c u = 1, 
and (c) c v = 1, c u = 1.99. The case (d) show the asymmetric dependence for A = 0.5 and c v = 1, 
c u = 1.7. 

mode with the energy inside the linear spectrum. Indeed, the values q(k min ) = and 
q(k max ) = Ti define the band edges of the propagation spectrum of the SH field, and such 
resonances take place when the SH field is generated. 



IV. CONCLUSIONS 



Being driven by the recent successful experimental demonstrations of discrete optical 
solitons with two frequency components mutually locked by a quadratic nonlinearity, in 
this paper we have studied two-color localized modes in arrays of channel waveguides. We 
have assumed that the tunable quadratic nonlinearity is introduced as a nonlinear defect by 
periodic poling of a single waveguide in the array, and we have analyzed the structure and 
stability of discrete localized modes created by mutual locking of two frequency components. 
We have shown that, depending on the value of the phase mismatch and the input power, 
such two-color defect modes can be realized in the same array in three different localized 
states, and we have studied also the resonant light scattering in the array with the defect 
waveguide drawing analogies with the Fano resonance. 
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